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Abstract 

We describe a high efficiency detector for measuring electron-positron pair transi- 
tions in nuclei. The device was built to be insensitive to gamma rays and to ac- 
commodate high overall event rates. The design was optimized for total pair kinetic 
energies up to about 7 MeV. 
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1 Introduction 



During stellar helium burning, the "triple alpha" (3a) reaction involving the 
excited 0"*" state at 7.65 MeV (the so-called "Hoyle state") in ^^C converts 
alpha particles to ^^C. Certain stellar processes are sensitive to the 3a reaction; 
current models of stellar evolution require that its rate must be known to 
within about 5% to describe these processes accurately [1]. The reaction rate 
is inversely proportional to the pair branch (F^/F), of the Hoyle state, the 
fraction of the time that, once formed, this state decays to the ground state 
by the emission of an electron-positron pair. This branch is presently known 
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Fig. 1. Schematic view of the detector. The letters correspond to those in Table 1. 
The Si detectors are not shown. 

with an accuracy of 9.2% [2,3,4] and is by far the largest uncertainty in 
determining the rate. 

The detector described here was designed to measure F^/F to within about 
5%, resulting in a determination of the 3a rate to about 6%. The experiment 
is simple in principle: the 0^ excited state in ^^C is excited by inelastic proton 
scattering at a bombarding energy of 10.4 MeV. The population of the state 
is tagged by the observation of scattered protons detected using two Si PIN 
diodes placed at 125° in the laboratory with respect to the beam direction. 
The choice of the scattering angle takes advantage of a resonance in the exci- 
tation function of the ^^C(p,p )^^C(7.65 MeV) reaction to maximize the count 
rate [5]. Positron-electron pairs from the decay of the 0^ state are detected 
in coincidence with protons using an array of plastic scintillators. The ratio 
of the proton-pair coincidence rate to the proton singles rate then gives the 
branching ratio, after correction for detection efficiency, random coincidences 
and gamma-ray contamination. The principal difficulty of this measurement 
is that F^/F is extremely small, about 6.7 x 10~^, and hence sensitive to even 
small backgrounds for the e"'"-e~ pairs. The dominant physical background is 
the much more probable cascade gamma decay of the Hoyle state which pos- 
sesses a branching ratio of F^/F = 4.1 x 10"'^. Excitation of the 2"*" state at 
4.44 MeV in ^^C and of states in target impurities, e.g. ^^Si, are other sources 
of background. Additional (random) backgrounds can be produced by inter- 
actions of the proton beam with collimators, elements of the detector, and the 



2 




Fig. 2. Photograph of the detector installed in the scattering chamber. Lead shields 
were placed in front of the beam entrance hole to reduce background due to the 
scattering of the beam on the collimators. 

beam stop. Further details on the experiment are given in References [1] and 
[6]. 

Section II describes the detector constructed for this purpose with special 
attention to features important for the reduction of backgrounds. In sections 
III and IV the results of GEANT4 simulations are discussed and compared 
with the observed performance of the detector. 



2 Detector properties and construction 

The conceptual design of the detector is shown in Figure 1. Figure 2 shows 
a photograph of the detector installed inside the scattering chamber. The 
active elements of the detector, consisting of a thin tube surrounded by a 
square block segmented into four quadrants, were constructed of Bicron BC404 
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Table 1 

Detector materials and dimensions-letters refer to Figure 1. All circular dimensions 
are diameters, with OD to mean "outer diameter" and ID, "inner diameter". 



Object 



Material 



Dimensions 



Scintillator block (a) 



Scintillator tube (b) 



Absorber tube (c) 



Light guides (d) 
Beam hole liner (e) 

Paint 
Quadrant PMTs (f ) 
Tube PMT (g) 
Si detectors 



Bicron BC404 



Bicron BC404 



Clear acrylic 



Bicron-800 
Tantalum 
F113 Epoxy with Ti02 
Electron tubes 9845B 
Hamamatsu R6231 
Si PIN diodes 



10x10x10 cm^, 
central hole 50 mm , 
four quadrants 
height 105 mm, 
40 mm ID, 46 mm OD, 
all surfaces polished 
height 105 mm, 
34 mm ID, 38 mm OD, 
all surfaces polished 

11 mm ID, 12 mm OD 
about 0.25 mm thick 
3.2 cm photocathode 
4.6 cm photocathode 
1 mm thick, 3 sub-sectors 



plastic scintillator. The energies deposited by the electrons and positrons in the 
quadrants of the scintillator block were read out using photomultiplier tubes 
(PMTs) attached to each quadrant with a light guide constructed of BC800 
plastic. The scintillator tube was coupled to a PMT using a flexible sihcone 
pad. By requiring coincident events in one or more quadrants and the central 
tube, gamma-ray backgrounds were suppressed by the low sensitivity of the 
thin scintillator tube to gamma rays. In order to shield the active elements 
of the scintillator array from protons scattered from the target, a 2 mm thick 
tube of inert plastic (acrylic) surrounded the target inside of the active tube as 
shown in Fig. 1. Holes were bored in this central tube to allow protons in the 
beam, or protons that were to be detected in the silicon detectors to emerge. 
Table 1 summarizes the materials used in the construction of the detector. 



The absorber and scintillator tubes, as well as the surrounding scintillator 
blocks, were mounted on an aluminum base plate. The surfaces of the scin- 
tillator quadrants were covered with white epoxy paint so as to provide light 
insulation between the quadrants, as well as to improve light collection. The 
segmentation of the square block into quadrants improves the rate capability 
of the device by limiting pulse pile up. 
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The beam entrance and exit holes were hned with tantalum tubes which pre- 
vent beam particles from entering the scintillators. The beam holes were also 
painted with white epoxy paint before insertion of the liners so as to prevent 
the absorption of scintillation light by the dark hner. 

The silicon detectors were 1 mm thick PIN diodes segmented into three sectors 
and are described in Reference [7]. The solid angle for protons is defined by 
two 4 mm diameter holes in the acrylic tube at 125° in the laboratory with 
respect to the beam direction. These holes extend through the scintillator tube 
with a diameter of 12.7 mm; the holes through the adjacent quadrants have 
a diameter of 12.7 mm, widening to 16 mm at the light guides to prevent the 
exiting particles from scattering back into the scintillators. The size of the 
apertures was chosen to keep the kinematic spread in proton energy compara- 
ble to the intrinsic Si detector resolution. The acceptance solid angle can be 
further restricted by tantalum apertures placed in front of the Si detectors, so 
as to reduce the kinematic spread. 



3 GEANT4 simulations 

We have carried out Monte Carlo simulations of the scintillator detectors using 
the code GEANT4 [8]. The aim of the simulations is to guide the design, as well 
as to understand the response of the scintillators to electrons, positrons, and 
gamma rays. This understanding is crucial to the determination of important 
operational characteristics such as pair efficiency, gamma-ray efficiency, and 
suppression. In practice, these properties are determined from comparisons 
between data and simulation for isolated, well-identified transitions. For the 
gamma-ray cascade of the -'^^0(02') state, the response to the 3.21 MeV gamma 
ray is determined from the Monte Carlo simulation. 

3.1 Elements of the GEANT4 simulation 

The GEANT4 simulation includes all the active elements of the plastic scintil- 
lator detectors. These elements are modeled with a realistic detector geometry. 
The central, inactive tube is also included in the simulation as particles and 
gamma rays deposit energy in this element as well. The simulation provides 
the energy deposited in each volume included in the experimental geometry. 
Elements not included are the light guides, aspects of light collection and light 
attenuation, and the response of the PMTs. Experimental resolution effects 
are treated later in the simulation process as described below. 

Electron-positron pairs from EO transitions were simulated using the well- 
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known energy and angular correlations given by the Born approximation [9]: 

^ = (e2 _ 1)1/2 X (^/2 _ X _ 

etc 

where e and e' are the total energies of the electron and positron. The angular 
correlation between the electron and positron is given by: 

d'^N dN 

— X (1 + a X cosB) (2) 



dedcos9 de 

where a is the EO anisotropy factor derived using the Born approximation (see 
Appendix C of Reference [10]), and is given by: 

^= 



The response of the detector to pairs from the decay of the ^^0(02"; 7.65 MeV) 
and ^^0(0^; 6.05 MeV) states was studied. The 0^ state in ^^O provides a 
convenient calibration signal as it has a 100% e"'"-e~ decay branch to the ^^O 
ground state. All particles were followed and their energy losses monitored 
until they either stopped, annihilated, or left the detector. 

The response of the detector to 4.44 and 3.21 MeV gamma rays (the energies 
produced in the gamma-ray cascade decay of the ^^C(02 ) level), gamma rays 
from radioactive ^°Co and ^^^Cs sources, and a variety of other transitions in 
^^Si and ^^O were simulated. In all cases, the photons were emitted from the 
target center and distributed isotropically throughout the entire solid angle. In 
the case of two-gamma-ray cascades, the photons were each emitted isotrop- 
ically and the angular correlation between them was ignored. For positrons, 
electrons, and gamma rays the low-energy processes of GEANT4 were used 
to extend the validity range of the coded electro-magnetic interactions and 
cover the full range of experimental energies (see Chapter 12 of Reference 
[8]). These processes include annihilation in flight for positrons. 



3.2 e^-e response and efficiency 



The design goals for the detector were that it be able to stop and measure the 
energies of the positrons and electrons produced from the ^^C(02 ) pair decay, 
while remaining relatively insensitive to the gamma-ray cascade decay of that 
state. Figure 3 illustrates the position within the detector, where positrons 
or electrons from EO transitions from the (a) ^^C(0^; 7.65 MeV) and (b) 
^^0(02"; 6.05 MeV) states stop. Each figure shows the results of 10^ e+-e~ 
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Fig. 3. The x and z coordinates (in mm) of the points where the electrons and 
positrons from the pair decay of the (a) ^^C(0^; 7.65 McV) and (b) ^^O(0j; 6.05 
MeV) states stop for 10^ generated pairs. The thick hnes show the contours of the 
four scintillator blocks surrounding the acrylic tube and the scintillator tube. For 
those particles which step outside the detector, that first step outside is shown. The 
positrons which annihilate in flight are not represented in the plots. 

pairs emitted at the target position in the center of the plot. Electrons and 
positrons stopping in the two central tubes as well as in the scintillator cube 
are clearly identified. On both figures the effects of the Ta liners for the beam 
entrance and exit can be seen (between x=±5 mm). 

Most of the electrons and positrons from the 7.65 MeV state in -"^^C which 
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Tube energy (MeVj Total energy (MeVj 



Fig. 4. (a) Simulated spectrum of energy deposited in the active scintillator tube 
by 10^ e"'"-e~ pairs from the decay of the ^^0(02"; 6.05 MeV) state, (b) Spectrum 
of total energy deposited in all active detector elements by pairs produced from the 
decay of the ^^O(Oj) state. 

enter the quadrants, stop near the center of the detector and away from the 
corners where light collection presumably deteriorates. Those positrons or elec- 
trons that exit the detector without depositing their full energy mainly do so 
through the central hole, the beam entrance and exit holes, the empty spaces 
between the acrylic and the scintillator tubes or the scintillator tube and the 
quadrants; some particles exit the detector after depositing energy in a quad- 
rant. 

Figures 4(a) and 4(b) show the energy deposited in the active scintillator 
tube and the total energy deposited in all active components of the detector 
(i.e., the scintillator tube and the four quadrants) for simulated ^^0(02 ) pair- 
decay events. The low-energy peak structure observed in Fig. 4(a) arises from 
events where one of the two particles never reaches the tube, either because 
it exited the device or was absorbed by the inert central tube. The higher 
energy peak corresponds to events where both particles deposit energy in 
the scintillator tube. In Figure 4(b), the spectrum peaks at a total deposited 
energy of about 4.18 MeV, showing that on average approximately 1 MeV of 
pair kinetic energy is lost by particles traversing the inert central tube. 

To determine the simulated detection efficiency for e+-e^ pairs, we make the 
following assumptions: for an event in which the energy deposited in the scin- 
tillator tube was greater than 0.16 MeV, and the energy deposited in at least 
one of the four scintillator quadrants was greater than 0.13 MeV, correspond- 
ing to the respective experimental energy thresholds, the pair was detected. 
Otherwise, the pair was not detected. The ratio between the number of de- 
tected pairs and total number of emitted pairs gives the simulated detection 
efficiency, and was ~85% for the decay of the ^^©(Oa") state, and ~90% for the 
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Fig. 5. Simulated energy-deposition spectra for 10'' generated events for the tube 
(left hand side) and all active detector elements (right hand side) for gamma rays 
with energy (a), (b):l. 17+1.33 MeV (^OCo); (c),(d): 1.78 MeV (^^Si); (e),(f): 4.44 
MeV (i2c); (g),(h): 6.92 MeV (^^O) and (i),(j): 7.12 MeV (^^O). The spike at 1.022 
MeV below the total energy corresponds to escape of both 511 keV annihilation 
gamma rays. 



^^0(02) state. Since this difference is small, observation of pairs from ^^0(02") 
provides a good baseline for an estimate of the experimental efficiency for pairs 
from the ^^C(02") excitation, which cannot be determined directly. While the 
experimental detection efficiency depends on other factors in addition to the 
total deposited energy, such as timing, detector thresholds, etc., the differ- 
ence between the efficiencies for the two pair transitions should be relatively 
insensitive to these detailed effects. 
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Fig. 6. Simulated and experimental detection efficiencies for single gamma-ray 
(28Si(1.78 MeV), i2c(4.44 MeV), 1^0(6.92 MeV), and i6o(7.12 MeV)) and cascade 
gamma-ray transitions (28Si(2.84 + 1.78 MeV), 28Si(3.20 + 1.78 MeV), 28Si(4.50 + 
1.78 MeV), and ^^81(6.02 + 1.78 MeV)). See the text for an explanation of the dis- 
crepancies between the simulated and the corresponding experimental efficiencies. 
For cascade gamma-ray transitions, the sum of the energies of the two gamma rays 
is shown on the x-axis. The error bars are purely statistical. Where not visible, the 
uncertainties are smaller than the size of the symbols. The statistical uncertainties 
for the simulations are negligible. 



3.3 Response to gamma rays 



Figure 5 shows the simulated total energy deposited in the scintillator tube and 
all active detector elements for several gamma-ray transitions as detailed in 
the caption. The sensitivity of the device to the 511 keV photons is very small 
and therefore an additional 1.022 MeV of the single gamma-ray transition 
energy is lost for many pair production events. The sharp peaks 1.022 MeV 
below the maximum deposited energy in the total energy spectra illustrate 
this effect. 

The efficiency for gamma rays estimated from these simulations was obtained 
in the same fashion as described above for c+-c^ pairs. The simulated efficien- 
cies for single gamma-ray and cascade gamma-ray transitions are shown in Fig- 
ure 6, where they are compared to their experimental values. The experimental 
gamma-ray coincidence efficiency was determined by comparing proton events 
obtained with, and without coincident signals in the scintillator-detector ar- 
ray. For a given total gamma-ray energy, the efficiency for detecting at least 
one of the two photons in a cascade is systematically larger than that for a 
single photon carrying the same energy. Figure 7 shows (a) a "singles" proton- 
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Fig. 7. Proton energy for one segment of one of the silicon detectors showing (a) 
the singles and (b) the coincidence spectra in the energy range 1.7-7.7 MeV for the 
Si02 target. The single gamma ray transitions for which the estimated experimental 
efficiencies are shown on Figure 6 are labeled in addition to the ^'^0(6.05 MeV) 
and the ^^0(6.13 MeV) transitions. All peaks have been identified; the majority are 
from Si. The statistics in (a) reflect a down-scale factor of 10. 

energy spectrum from pH-Si02 collisions obtained without the requirement of 
a scintillator coincidence, and (b) the corresponding spectrum obtained with 
silicon-detector signals in coincidence with pulses in the scintillator detectors, 
corrected for random coincidences. Several proton-energy peaks are observed 
corresponding to different excitations in ^^O and ^^Si, as well as in ^^C that 
was present as an impurity in the target. Clearly visible in Figure 7(b) is the 
relative enhancement of the peak of the pair-emitting ^^0(62") excitation at 
Ep = 3.086 MeV compared to excitations that decay by gamma-ray emission. 
The experimental efficiency for e+-e~ pairs is discussed below. 
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Fig. 8. Two-point tube calibration using (a) the peak at 0.5 MeV for ^^0(4.44 MeV) 
and (b) the one at 1.2 MeV for ^60(6.05 MeV). The thick dark hues show a fit to 
the data using an empirically determined detector resolution function for the tube 
(see text). 

4 Comparison Between Experimental Results and Simulations 

4-1 Energy calibrations 

Figure 8 illustrates the method used to calibrate the tube energy. The peaks in 
the tube spectra near 0.5 MeV (^^C{AAA MeV)) and 1.2 MeV (^^0(6.05 MeV)) 
are first assigned a reasonable channel number based on the corresponding 
uncalibrated experimental spectra. The calibration constants are then refined 
recursively by fitting the simulated spectra to data. Detector thresholds and 
an approximate detector resolution function are applied to the simulations 
for each detector component based on experimental data. Once the energy 
calibration has been determined, the best fit parameters for the resolution 
function used to smear the simulated data are obtained by fitting the tube 
spectrum for ^^0(6.05 MeV). They are then fixed to fit the tube spectra for 
many gamma rays (we illustrate the case of ^^C(4.44 MeV) on Figure 8(a)). 
The same method is used to determine the best fit parameters for the quad- 
rants. This gives the following detector resolution functions: a Gaussian with 
an energy dependent width of a{MeV) = 0.00097 + 0.16 x VE + 0.0010 x E 
for the tube spectra and a{MeV) = 0.0045 + 0.088 x + 0.00098 x E for 
the quadrant spectra, where E is the deposited energy. The widths of the res- 
olution functions include three terms: an energy-independent part to account 
for the baseline noise, a ^/E term to account for the statistical nature of light 
collection, and a term proportional to E to account for the variation of light 
collection efficiency within the detector. Figures 9(b) and 9(d) show the same 
two-point calibration technique applied to one quadrant using the ^^C(4.44 
MeV) and the ^°Co(l. 17+1.33 MeV) transitions. The high-energy drop-off on 
the quadrant spectra occurs below the Compton edge for gamma rays from 
^^C and ^^Si due to the requirement of energy being also deposited in the 
scintillator tube. The simulations for Si (1.78 MeV) are convoluted using the 
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Fig. 9. Left: unsmeared simulated 1-quadrant gamma-ray spectra. Right: the cor- 
responding experimental spectra for gamma rays of energy (a),(b):1.17+1.33 MeV 
(60Co), (c),(d):4.44 MeV (^^C), (e),(f):1.78 MeV (^^Si). A threshold of 0.16 MeV for 
the tube and of 0.13 MeV in at least one of the quadrants has been applied to the 
simulated spectra. The thick dark lines show a fit to the data using an empirically 
determined detector resolution function for the quadrants (see text). 

same resolution function and normalized to data as shown on Figure 9(f). The 
results are in good agreement with the data. 



Figure 9 also shows a comparison between the simulated (unsmeared) and 
the random-subtracted data spectra for various gamma-ray transitions after 
energy calibration has been performed. Figures 10(a) and 10(b) show the un- 
smeared simulated spectra for energy deposited by e+-e~ pairs for the decay 
of ^^0(6.05 MeV) in one quadrant and all four quadrants, respectively. Reso- 
lution smearing is not included in Figures 9(a), 9(c), 9(e), and Figures 10(a) 
and 10(b); resolution smearing is, however, performed before normalization 
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Fig. 10. Spectra for (a),(b) the unsmeared simulations and (c),(d) the experiment 
for energy deposited by electron-positron pairs from the ^^O(O^) in 1 quadrant 
((a),(c)) and in all 4 quadrants ((b),(d)). In (c) and (d), the thick dark line shows 
a fit to the experimental spectra of the corresponding simulated spectra smeared 
with the appropriate detector rcsohition function (see text). 

to data. Figures 10(c) and 10(d) show the corresponding experimental energy 
spectra normalized to the simulations. The same detector resolution function 
as the one used for the gamma-ray quadrant spectra has been applied to the 
simulated spectra before normalization; the results are in good agreement with 
the data. 



4-2 e^-e detection efficiency and gamma-ray suppression 



The experimental efficiency of the scintillator-detector array for e+-e~ pairs 
from the ^^0(02^) level was obtained in a manner similar to that described 
above for gamma rays. Fig. 11 shows proton-energy spectra in the region 
of the ^^0(02") peak for (a) singles proton events and (b) random-subtracted 
coincidence events. A complication in the measurement arises from the nearby 
3- excitation in ^^O at Ep^ 3.015 MeV and a 0+ excitation in ^^Si at Ep= 2.99 
MeV, which are not completely resolved from the ^^0(02") peak. To extract 
the yield for the ^^O(O^), the data were fit with three overlapping Gaussians. 
The curves in Fig. 11 show the resulting fits to the proton-energy data. In 
Fig. 11(a) the sum of the 3" (^^O) and 0+ (^^Si) peaks is dominant. With the 
application of the scintillator-detector coincidence the ^^0(02") peak is strongly 
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Fig. 11. Proton-energy spectra from one segment of one of the silicon detectors for 
(a) singles and (b) coincidence events for the Si02 target. The singles spectrum 
has a down-scale factor of 10. The individual components of the fit include a linear 
background and three Gaussian peaks corresponding to ^^Si(6.69 McV) at Ep=2.99 
MeV, ^^0(6.13 McV) at Ep=3.015 MeV, and ^^0(6.05 MeV) at Ep=3.086 MeV. 

enhanced relative to those corresponding to excitations that decay by gamma 
rays. The experimental value of the pair efficiency is 65.4 ± 1.8 %, compared 
to the simulated value of 85%. The systematic over-prediction of the detector 
efficiency is similar to that observed for various gamma-ray transitions.The 
difference between the observed and simulated efficiencies is attributed to 
effects not treated in the simulation, such as light collection in the scintillators 
and light guides and a realistic treatment of electronic thresholds. While these 
effects are not fully addressed in the simulation, the experimental efficiency 
may be calibrated by the ^^0(02 ) pair measurement, and extrapolated to the 
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nearby ^^C(O^) state with guidance from the Monte Carlo calculation. 

5 Conclusion 

The pair detector described here was used for an experiment to better deter- 
mine the rate of the triple alpha experiment by improving the error bars on 
the internal pair emission branch of the 0^ excited state in ^^C. The detector 
performance was simulated using the GEANT4 simulation framework. The 
simulations helped determine the best detector geometry at the design stage 
and were crucial for the detector calibration and the determination of impor- 
tant detector efficiencies. Monte Carlo simulations combined with empirically 
determined detector resolution functions yielded good agreement with the ex- 
perimental data. 
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